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ABSTRACT 

The parameters of the GRB 000301c afterglow are determined within the usual framework of synchrotron 
emission from relativistic ejecta, through fits to the available radio-to-optical data. It is found that 1) the jet 
energy after the GRB phase is ^ 3 x 10^" erg, 2) the initial opening angle of the jet is ~ 13°, 3) the medium 
that decelerates the afterglow has a density ^ 25 cm~^, 4) the power-law distribution of the shock-energized 
electrons is hard, with an index around 1.5, and 5) the cooling frequency was located below the optical domain. 
Furthermore we find that the collimation of ejecta alone cannot explain the sharp and large magnitude break 
observed in the i?-band emission of this afterglow at few days, and that a high energy break in the electron 
distribution, corresponding to an electron energy close to equipartition, together with the lateral spreading of the 
jet accommodate better this break. Microlensing by a star in an intervening galaxy produces only a flattening of 
the afterglow emission and cannot explain the mild brightening exhibited by the optical emission of 000301c at 
~ 4 days. 

Subject headings: gamma-rays: bursts - ISM: jets and outflows - methods: numerical - radiation mechanisms: 
non-thermal - shock waves 



1. INTRODUCTION 

The afterglow of the GRB 000301c is one of the best ob- 
served to date, its emission having been monitored over 1.5 
decades in time at radio frequencies (Berger et al. 2000) and 
in the optical (Bhargavi & Cowsik 2000, Masetti et al. 2000, 
Sagar et al. 2000, Jensen et al. 2001, Rhoads & Fruchter 2001). 
Furthermore, 000301c is one of three afterglows whose opti- 
cal emission decay exhibited a strong steepening (break) at few 
days. The steepening follows a flattening, or mild brightening, 
in the optical light-curve. The break in the optical emission was 
attributed to a jet geometry by Berger et al. (2000), while the 
relatively short-Uved brightening at 3-5 days was interpreted by 
Garnavich, Loeb & Stanek (2000) as evidence for gravitational 
microlensing due to a star in a intervening galaxy. Rhoads & 
Fruchter (2001) point out that the sharpness of the optical break 
poses difficulties for a jet interpretation. 

In this work we model the emission of 000301c in the frame- 
work of relativistic ejecta energizing the swept-up gas sur- 
rounding the GRB source (Meszaros & Rees 1997). Our aims 
are to assess the ability of the jet model to accommodate the be- 
haviors exhibited by the radio-to-optical emission of 000301c 
and to determine its basic physical parameters (jet energy, ini- 
tial opening, surrounding medium density, etc). Furthermore, 
by including microlensing in our model, we also test its ability 
to reproduce the optical brightening seen around 3.5 days. 

2. MODEL DESCRIPTION 

2. 1 . Afterglow Dynamics and Synchrotron Emission 

The jet dynamics is given by a set of differential equations 
(Panaitescu & Kumar 2001) for the Lorentz factor F, mass, 
and aperture 6 of the jet as functions of radius r. The jet lat- 
eral spreading is negligible in the early phase of the expan- 
sion. In the adiabatic limit, the jet Lorentz factor evolves as 
r cx r~^/^ oc t~^^^ for a uniform external medium, where t 
is the observer time, and as F cx r~^/^ cx t^^/* for a medium 
whose density decreases as n{r) cx r^^ (e.g. Meszaros, Rees & 
Wijers 1999), i.e. for a wind ejected at constant speed before the 
release of the relativistic jet. The lateral spreading increases the 



jet sweeping area and alters significantly its dynamics when F 
drops below (Rhoads 1999), the jet deceleration becoming 
exponential, F cx e"*^"" cx t~^^^, as long as the jet is sufficiently 
relativistic (F > 10). 

We calculate numerically the jet dynamics, taking into ac- 
count the energy loss through emission of radiation. The jet is 
assumed to have sharp boundaries and to be uniform, with no 
angular gradients (of energy or mass) within its aperture. From 
the dynamics of the jet, one can calculate the properties of the 
shocked fluid (density, magnetic field, typical electron random 
Lorentz factor), which are necessary for the computation of the 
synchrotron emission. 

We assume that the shock injects in the downstream region 
relativistic electrons with a power-law distribution, and approx- 
imate the distribution of electrons resulting from injection and 
cooling as a dual power-law (Sari, Piran & Narayan 1998), with 
a minimum and a break set by the injected (7^) and cooling (7c) 
electron Lorentz factors. This approximation allows some an- 
alytical results for the self-absorption frequency, Compton pa- 
rameter, radiative losses, and peak synchrotron flux to be im- 
plemented in the code, avoiding thus time-consuming integra- 
tions over the electron distribution. We have tested the accuracy 
of this approximation and find that it is very good if radiative 
losses are small. For significant radiative losses (above 50%) 
it may mis-estimate the electron distribution and 7c by a fac- 
tor ^ 3, thus the model fluxes in some frequency ranges at for 
certain durations may be off by a similar factor. 

In the calculation of the afterglow emission we allow for the 
existence of a high energy break at 7* in the electron distribu- 
tion, at which the shock-acceleration becomes inefficient. This 
break could be due to the escape of particles from the accel- 
eration region and to radiative cooling during acceleration, and 
becomes important if the emission from 7, -electrons crosses an 
observational band. In this case the resulting afterglow emis- 
sion fall-off depends on the exact shape of the electron distri- 
bution at and above 7,^. For simplicity we shall assume that at 
7, there is a sharp transition to a steeper power-law. 

The afterglow spectrum has three breaks: the synchrotron 
characteristic frequencies i^i and i^c at which the 7^- and 7c- 
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electrons radiate, and the self-absorption frequency Va- Be- 
tween these breaks the spectrum is a power-law of an index de- 
pendent on the ordering of the breaks (e.g. Sari et al. 1998). The 
steepening of the electron distribution at 7* adds another break, 
at the frequency v^, at which the 7* -electrons radiate. As the jet 
is decelerated, these four frequencies evolve, the passage of one 
of them through an observing band yielding a chromatic steep- 
ening of the afterglow emission fall-off. Furthermore an achro- 
matic steepening of the afterglow decay is expected around the 
time when, due to the decrease of the relativistic beaming of the 
afterglow emission, the jet edge becomes visible to the observer 
(Rhoads 1999). 

In our treatment, the afterglow model has 8 parameters. 
Three of them give the jet dynamics: 

1) Eq, the jet energy at the end of the GRB phase 

2) Oq, the initial jet aperture (half-angle) 

3) n, the particle density of the uniform external medium or 
the parameter A = (l/47r) (M/nipv) for a medium with a den- 
sity profile n(r) = Ar~^, resulting from the ejection by the 
GRB progenitor of a wind with mass loss rate M and constant 
speed V. The other five are parameters for the microphysics of 
shock acceleration and generation of magnetic fields, and de- 
termine the spectrum of the emergent synchrotron emission: 

4) £b < 1/2, the ratio of the magnetic energy density to the 
energy density of the freshly shocked fluid. 

5) £e 1, the ratio of the energy density of newly acceler- 
ated electrons, if they all had the minimum Lorentz factor 7j, to 
the energy density of the post-shock gas. This parameter quan- 
tifies 7i. 

6) p, the index of the power-law distribution ^(7) oc 7"^ 
(for 7i < 7) of the injected electrons. 

7) e < 1/2 (greater than Eg), the fractional internal energy 
stored in electrons. This parameter determines the high-energy 
break 7* of the electron distribution for p < 2 and any Sg < e 
or for p > 2 and e < [{p - l)/ip- 2)]£e- 

8) Sp, the steepening of electron distribution index above 7*. 
A quantitative description of these five parameters can be found 
in Panaitescu & Kumar (2001). 

We note that, if the jet evolution is quasi-adiabatic, then its 
dynamics and the afterglow emission are independent on Fq. 
If radiative losses are significant, the afterglow emission de- 
pends on Fo because the magnitude and evolution of the radia- 
tive losses is Fo-dependent. However, this dependence is rather 
weak, so that the data cannot be used to obtain meaningful con- 
straints on the initial jet Lorentz factor. For this reason we shall 
keep it fixed: Fq = 500. 

In the calculation of the radiative losses and cooling fre- 
quency we take into account synchrotron self-absorption and 
inverse Compton scatterings (which may be dominant if the 
magnetic field is sufficiently weak). The received synchrotron 
flux is calculated by integrating it over the jet evolution, tak- 
ing into account the effect of the source relativistic motion and 
spherical curvature on the photon arrival time and energy, but 
ignoring the light shell-crossing time. The observer is assumed 
to lie on the jet axis. 

2.2. Interstellar Scintillation 

Fluctuations on a short timescale of the radio emission of 
GRB afterglows, caused by inhomogeneities in the Galactic 
interstellar medium have been predicted by Goodman (1997). 
As shown by Berger et al. (2000), the weU-sampled 8.46 GHz 
emission of 000301c exhibits such fluctuations. We include the 



effect of the interstellar scintillation (ISS) by adding in quadra- 
ture to the reported errors the expected amplitude of the model 
flux fluctuations due to ISS. The calculation of the modulation 
amplitude follows the treatment given by Walker (1998), and is 
summarized below. 

At a frequency v below a certain critical value fo, scin- 
tillation occurs in the strong regime. Diffractive scintiUa- 
tion (DS) is caused by the interference of Ught coming from 
many patches of coherence within the "refractive" disk of size 
6r = Opoij^o/i^y^^^, where 6*^0 is the angular size of the first 
Fresnel zone at the critical frequency vo- DS has amphtude 

m'f^^ = 1 if the source is point-like (i.e. its size 9s is smaller 
than 0d = OfoW/i^o]^^^) and if the observing window is nar- 
rower than the decorrelation bandwidth Au = i/{v/voy^/^. 
For an "extended" source (Od < 9s) the modulation index is 

reduced to mll^^ — Od/Og. 

Refractive scintiUation (RS) is due to inhomogeneities in the 
scattering screen on the refractive disk scale. In the strong 

scattering regime the RS amplitude is mr^*'' = (iz/i/g)^''''^'^ 
if 9s < 9r, while for 9r < 9s the modulation index is 
mi''^^ = mi^^\9r/9sy^^. At frequencies above uq, RS oc- 
curs in the weak scattering regime and have amplitude 
[vQlvf^l^"^ if 61,, < 9f, where 9f = 9Fn{i^n/vY'^ is the size 
of the Fresnel zone. For 9f < 9s the RS amplitude is reduced 

to m^^'^ =m^^'\9F/9sy/^ 

The amplitude of the radio ISS depends therefore on two 
properties of the interstellar medium, Vq and 9fo7 which de- 
pend on the scattering measure and distance to the scattering 
screen. The maps of Walker (1998) give 9fo ~ 5.0 /zas and 
uq ~ 4.4 GHz in the direction of GRB 000301c, thus the mea- 
surements at the two lowest frequencies (4.86 GHz and 8.45 
GHz) are close to the critical frequency. The above equations 
for the ISS amplitude are asymptotic results, strictly valid at 
frequencies far from fQ, and slightly overestimate the modula- 
tion amplitude atv uq. For this reason we set an upper limit 
of 50% to the modulation index calculated with the asymptotic 
results above. 

The source size 9s, on which the ISS modulation may de- 
pend, can be easily calculated analyticaUy at times before lat- 
eral expansion becomes important. During this phase the jet 
aperture 9 exceeds F~^, thus the observer receives photons 
mostly from a spherical cap of opening 9 = around the 
jet axis. Therefore the size of the source seen by the observer 
is 9s = r/r. For a uniform medium (F cx r^^/^) it can be 
shown that a photon emitted by the fluid moving off-axis at an 
angle 9 arrives at observer 5 times later than of a photon travel- 
ing along the jet axis. For a wind-like medium (F oc r~^/'^) the 
above factor is 3. Applying these correction factors to the ex- 
pressions for the fireball radius given by Panaitescu & Kumar 
(2000), which contain the arrival time of the photons emitted 
along the jet axis, we obtain 

9s = 0.89 D^^^{1 + z)-^'\£o,5z/naY'^ Mas (1) 
for a homogeneous medium and 

9s = 0.60 D^]^^{1 + ^)-^/^(5o,53/A*)'/^ iT Mas , (2) 

for a wind-like medium, where f = 2 i?o/(l — cos 9o) is the 
jet isotropic equivalent energy and ^* = A/(3.0x 10^^ cm~^) 
is A normalized to the value corresponding to the ejection of 
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10"^ Moyr^i at 10^ km s^^). In equations ^ and (|) Da is 
the angular distance for redshift z, is the observer time mea- 
sured in days, and the usual scaling Xn — 10^"X was used. 

After r drops below 9^^, which is roughly the onset time 
of the lateral spreading dominated phase, the source size is 
9s = rsinO. This quantity and the arrival time of photons 
emitted from the jet edge are calculated numerically. 

2.3. Gravitational Microlensing 

The decay of the near infrared (NIR) and optical emission of 
000301c exhibits a modest brightening at 3-5 days, which ap- 
pears to be achromatic in the optical (Berger et al. 2000, Bhar- 
gavi & Cowsik 2000), although there is evidence for R — K 
color evolution during it (Rhoads & Fruchter 2001). Garnavich 
et al. (2000b) have proposed gravitational microlensing as an 
explanation for this feature, though other possibilities - inho- 
mogeneities in the external medium, delayed energy injection, 
departures of the electron distribution from a power-law - can- 
not be ruled out given the lack of well sampled, simultaneous 
measurements in other frequency domains and of a detailed 
analysis of the chromaticity of the brightening in all these sce- 
narios. 

We include microlensing in our modeling to assess its abil- 
ity to explain the brightening of the 000301c afterglow. Unlike 
the other scenarios mentioned above, lensing is "external" to 
the afterglow modeling, as it does not modify the jet dynamics. 
Better yet, it introduces only two free parameters: 

1) d^s^ the apparent separation between the lens and the 
center of the GRB remnant and 

2) 9e, the angular size of the Einstein disk, which depends 
on the lens mass and distance to the observer Thus the total 
number of parameters of the afterglow model with lensing is 
10. 

As pointed out by Loeb & Perna (1998), a time-varying 
magnification of afterglows results when the expanding source 
(eqs. [|l]| and [^) crosses the Einstein disk. The effect is en- 
hanced by that the source has a non-uniform surface brightness 
distribution, with most of the optical emission coming from the 
outer part of the visible disk. These two factors are taken into 
account in our calculations by integrating over the jet surface 
the flux received from each infinitesimal patch multiplied by 
the magnification factor 



u\J V? + 4 



(3) 



where u is the angular separation (on the sky) between the patch 
and the lens, measured in units qIBe- 

We note some important differences between the model used 
by Garnavich et al. (2000b) and ours. We calculate the after- 
glow emission from the jet dynamics instead of approximat- 
ing it with a smoothed broken power-law. The source surface 
brightness distribution is implicitly taken into account by our 
integration of the afterglow emission. Thus we do not approxi- 
mate the afterglow image as a ring, with zero brightness inside 
the ring. Garnavich et al. (2000b) used a free normalization 
factor for the afterglow emission in each of the 7 observational 
bands included in their fit. This allows a rather excessive free- 
dom to their model, given that the afterglow spectrum in NIR- 
optical is a power-law, thus the normalizing factor at one fre- 
quency determines the factors for all other bands. All these dif- 
ferences point to that our model has significantly less freedom 
than that of Garnavich et al. (2000b). 



3. SPECTRAL/TEMPORAL FEATURES OF THE 00030 IC AFTERGLOW 

The afterglow data which we model are taken from Bhargavi 
& Cowsik (2000), Berger et al. (2000), Masetti et al. (2000), 
Sagar et al. (2000), Jensen et al. (2001), and Rhoads & 
Fruchter (2001). Some of these articles list measurements re- 
ported in GCN Circulars by Bernabei et al. (2000), Gal- Yam 
et al. (2000), Garnavich et al. (2000a), Halpern et al. (2000), 
Kobayashi et al. (2000), Mujica et al. (2000), and Veillet & 
Boer (2000). We have excluded the earliest i?-band magnitude 
reported by Sagar et al. (2000), which is inconsistent at the 
6cr level with an almost simultaneous measurement by Bhar- 
gavi & Cowsik (2000), and the UV flux reported by Smette 
et al. (2001), which could be affected by intergalactic Lya ab- 
sorption. We have averaged measurements separated by less 
than 0.01 day, whether they were independent observations or 
the same data analyzed and reported by different workers. The 
final data set contains 140 points. 

The NIR and optical magnitudes have been converted to 
fluxes using the photometric zero points published by Campins, 
Rieke, & Lebofsky (1985), and Fukugita, Shimasaku & 
Ichikawa (1995). To account for magnitude-to-flux conversion 
uncertainty we add 5% in quadrature to the reported uncertain- 
ties, noting that most of these uncertainties are equal to or ex- 
ceed 5%. The NIR and optical fluxes have corrected for Galac- 
tic extinction with E(B — V) — 0.053 (e.g. Rhoads & Fruchter 
2001), but no correction was made for a possible reddening in 
the afterglow's host galaxy. 

Throughout this work we denote by a and /3 the indices of 
the power-law afterglow light-curve and spectrum: 



(4) 



Over various ranges in frequency and times a and /3 are con- 
stants. 

3.1. Spectral Slopes 

For 000301c, the measurements made around t = 4.3 day 
span the widest frequency range. At the same time the optical 
emission exhibits a brightening. Figure 1 shows the afterglow 
spectrum at < ~ 5.0 day, when the optical bump seems to have 
faded, obtained from the data closest to this time, using inter- 
polations or extrapolations where necessary. 

The radio spectrum of 000301c has a slope (3r — —0.75 ± 
0.15, between the values expected for optically thin {f3r = 
— 1/3) and thick (/3r — —2) emission, thus the absorption fre- 
quency i/a must be around 10 GHz at 5 days. 

Using the 250 GHz measurements reported by Berger et 
al. (2000) and the /-C-band data of Rhoads & Fruchter (2001), 
we find that the millimeter-NIR spectral slope (3rnK is 



==0.67 ±0.04 at t = 5d 



(5) 



Rhoads & Fruchter (2001) found that, after correction for 
Galactic extinction only, the spectral slope between the K- and 
i?-bands is (see their Table 3) 



/3. 



KR 



0.72 ±0.04 t = 2.1d 

0.91 ±0.04 i = 3.0d 

0.99 ±0.15 i = 5.0d 

0.69 ±0.10 i = 7.6d 



(6) 



while at an earlier time, < = 1.8 day, the R — K color presented 
in their figure 2 gives I3kr — 0.47 ± 0.11. Rhoads & Fruchter 
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(2001) find that, in order to reconcile the data with no evolu- 
tion of Pkr.^ a systematic error of 0.08 mag has to be added 
in quadrature to the K and i?-band fluxes, and conclude that at 
least some of K — R color evolution is real. 

From spectroscopic observations, Feng et al. (2000) found 
that the optical spectral slope is /3o ~ 1.1 at t = 2 day, while 
Jensen et al. (2001) arrived at 



/3o = 



1.15 ±0.26 t = 4.0d 
1.43 ±0.28 t = 5.0d 



Equations (^ and (0) lead to 



0.23 ±0.26 t 
0.44 ±0.32 t 



4d 
5d 



(7) 



(8) 



The above results provide only weak evidence for a curved 
(even after dereddening for Galactic extinction), softening af- 
terglow spectrum. Jensen et al. (2001) and Rhoads & Fruchter 
(2001) attribute this departure from the power-law spectrum ex- 
pected in the fireball model to extinction in the host galaxy. We 
note that softening of the NlR-optical spectrum, if real, sug- 
gests an intrinsic feature of the afterglow, as caused by a spec- 
tral break whose frequency decreases, passing through the opti- 
cal domain. Given that the fastest possible decrease of the cool- 
ing frequency is oc (see Panaitescu & Kumar 2001) 
and that the steepening A/3 = 1/2 of the spectrum across I'c 
is smooth, it seems rather unlikely that the passage of would 
produce a detectable afterglow softening over a short timescale. 
This suggests that the evolving break is the j/* associated with 
the high end of the electron distribution (7*), which evolves as 
fast as i^i and can be sharper than the i^c -break. 

It is also possible that the chromatic evolution of the after- 
glow between 2 and 5 days is caused by the same mechanism 
that yields the optical brightening seen at ^ 3.5 days. Most 
likely this excludes microlensing as the brightening mecha- 
nism, as lensing should not produce a chromatic afterglow be- 
havior over less than a decade in frequency (note that, in gen- 
eral, microlensing of afterglows leaves chromatic signatures be- 
cause the surface brightness profile of the source is wavelength- 
dependent). 

3.2. Light-Curve Decay Indices 

The 8.46 GHz radio emission of the 000301c afterglow, 
shown in Figure 2, exhibits a power-law decline with 



= 1.0 ±0.2 at t>30d 



(9) 



A steeper and highly uncertain decline, between t^^ and i^^^, 
is seen in the 4.86 GHz emission. Due to the curvature of the 
emitting surface, relativistic jets cannot produce significantly 
different decay indices at frequencies that are a factor of only 
two apart, so we presume that the light-curve indices at 4.86 
GHz and 8.46 GHz are in fact the same and that the discrep- 
ancy can be resolved by interstellar scintillation. 

The millimeter decay index at 250 GHz is not well deter- 
mined by the data, nevertheless the 2<t upper limits shown in 
Figure 2 indicate that > 0.7. This shows that after few 
days the injection frequency Vi is below ^ 100 GHz, unless the 
external medium is wind-like and the ordering of frequencies is 
I'a < 250 GHz < I'c < Vi, in which case it is analytically ex- 
pected that a„i = 2/3 (Panaitescu & Kumar 2000), marginally 
consistent with the observations. 



The NIR and optical emission of 000301c exhibited a break 
at few days. Using a smoothed dual power-law fit, Rhoads & 
Fruchter (2001) found a steepening from an asymptotic aKi = 
0.09 at early times to aK2 = 2.26 at late times, in the K band, 
and from am — 0.69 to apa — 2.77 in the R band. Using 
a broken power-law approximation, Jensen et al. (2001) found 
aoi = 0.73±0.27andao2 = 2.67± 0.51. At late times, when 
the effect of microlensing should be negligible, Garnavich et 
al. (2000b) found ao2 — 2.9, close to the value reported by 
Sagar et al. (2000): ao2 = 2.97 ± 0.04. With a more complex 
fitting function and without the last two R magnitudes reported 
by Fruchter et al. (2000), which constrain the late time decay 
index, Bhai-gavi & Cowsik (2000) found = 0.70 ± 0.07 
and apa — 2.44 ± 0.29. Thus the asymptotic indices of the 
optical decay of 000301c are 



ttfli = 0.70 ±0.07 , afl2 > 2.7 



(10) 



implying that the decline of optical emission of this afterglow 
steepened by 

Aafl = afl,2 - a_Ri > 1.9 (11) 

over about one decade in observer time. 

For a homogeneous external medium, the afterglow decay in- 
dex is expected to increase by 3/4 when the jet edge becomes 
visible (i.e. the jet Lorentz factor drops below the inverse of 
its aperture, 9), with an extra steepening of at most p/A arising 
from the jet dynamics. It is thus hard to see how collimation of 
ejecta could account alone for the magnitude Aai? of the opti- 
cal break of 000301c, unless the afterglow decay before 3 days 
is significantly flattened by the brightening mechanism we see 
at 3-5 days, so that the observed decay index a^x is smaller 
than that the afterglow would have if this mechanism were not 
present. Given that it takes more than two decades in time to a 
jet expanding into an r^^ medium to yield about 90% of the full 
break magnitude (Kumar & Panaitescu 2000), the /S.aji above 
rules out a model with a wind-like medium and an afterglow 
steepening due only to collimated ejecta. 

4. MODELS FOR THE 000301C AFTERGLOW 

We start with a model where the high energy spectral break 
z/* is well above the optical domain and the light-curve steep- 
ening is due solely to the collimation of ejecta. In this case the 
external medium must be homogeneous, so that the light-curve 
steepening at the time tj, when the jet edge becomes visible, is 
as fast as possible. 

If the cooling frequency were below the NlR-optical domain, 
then (ixR = p/2, thus the observed slope at 5 days (eq. |Q|), 
after the optical brightening, would require p = 1.98 ± 0.30. 
Then the decay index of the optical emission at < > tj would 
approach asymptotically 032 ^ P S5 2.3, and would be shal- 



lower than observed (eq. |10|). Therefore the cooling frequency 
should be above the NlR-optical domain. 

At few days, the injection frequency Vi must be below the 
millim eter range, to explain the decay of the 250 GHz emission 
(§ [3.2| , Figure 2). Then (3mK — (p ~ l)/2 and the observations 
at 5 days (eq. [|l|) lead to p = 2.34 ± 0.08, which is still below 
C(R2- Therefore it may be difficult to reconcile within the jet 
model the spectrum and the late time optical decay of 000301c. 

The best fit obtained with a jet model having the cooling fre- 
quency above optical and a homogeneous external medium is 
shown in Figure 2. Its = 680 for 132 degrees of free- 
dom (df) is large in part because microlensing does not produce 
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the brightening seen around 3.5 days, for reasons discussed be- 
low, the 28 NIR/optical measurements between 3.0 and 4.3 days 
yielding Ax^ — 210. To test the ability of only the jet model to 
account for the general features of 000301c, we eliminate the 
3.0-4.3 days data and microlensing from the modeling. The 
new best fit has ~ 480 for 98 df (see Figure 2). 

One reason for which the new is still large is that a rel- 
ativistic jet yields the same decay after tj at all frequen- 
cies above vi. However, in 000301, the late time decay indices 
of the 8.46 GHz and i?-band light-curves differ significantly. 
For the jet models shown in Figure 2, the 14 measurements 
made at 8.46 GHz after 30 days yield A^^ — 68. Berger et 
al. (2000) explain this discrepancy (see their Figure 2) between 
the jet model and the radio emission of 000301c as the result 
of interstellar scintillation. Even if ISS were strong enough to 
accommodate the observations, the radio data should have been 
randomly scattered around the model light-curve, both having 
the same trend, which is unlike the behaviors illustrated in Fig- 
ure 2. 

The results of Figure 2 point to yet another reason for the 
large we obtain: the optical steepening of 000301c is too fast 
and too large to be explained with a j et b reak, as was expected 
based on analytical considerations (§3.2). For this reason we 
turn now to a "modified jet model", where the passage through 
the optical of the i/* break associated with the steepening of the 
electron distribution at 7* yields a light-curve steepening. As 
noted in § 3. 1 , the passage of a spectral break is suggested by the 
softening of the NIR-optical spectrum (eq. [||]). The time and 
magnitude of the light-curve break caused by the j/* -crossing 
depend on the location of 7* and on the steepening of the elec- 
tro n dis tribution at 7*, i.e. the parameters e and bp introduced 



in §2.1 



First note that, within the jet model, the quasi-constant 8.46 
GHz emission prior to 30 days has two possible explanations: 

1) the ejecta are coUimated, the jet edge becoming visible 
at few days, so that cx t^^l'^ before the passage of vi and 
Fy cx t^'P after that. These analytical results (Rhoads 1999) 
are not very accurate for moderately relativistic jets (F < 10), 
where we find numerically smaller decay indices. Nevertheless 
a hard electron distribution, with p < 1.5, is required by the ob- 
served ttr ~ 1. Since the observed PmK (eq- 0) is larger than 
(p — l)/2, the cooling break must be below NIR. The external 
medium can be either homogeneous or wind-like. 

2) the external medium is a wind and the ejecta spherical, the 
radio emission evolving as F^, w const before Vi — v and as 
F^j cx t^'y^'P-^)!^ thereafter. The observed (eq. [^) requires 
thatp = 1.7±G.3. As inthefirstcase, > (p—l)/2 implies 
that V(. is below NIR. Given that the jet break is very smooth if 
the external medium is wind-like, this case is not readily distin- 
guishable from a jet interacting with a wind. 

Figure 3 shows the best fits we obtain with a modified jet in- 
teracting with a homogeneous medium, with and without lens- 
ing and the 3.0^.3 day data. In the latter case, the resulting 
fit has = 120 for 96 df, being marginally acceptable. Note 
that, in either case, the modified jet yields a significantly better 
fit than the "traditional" jet model with p > 2 and optical steep- 
ening caused by coUimation. The best fits obtained with a mod- 
ified jet and a wind-like medium are shown in Figure 4. Without 
microlensing and the 3.0^.3 days data, it has x^ = 140 for 96 
df, being unacceptable. The light-curves shown in Figures 3 
and 4 prove that, for either type of external medium, microlens- 
ing does not yield an afterglow brightening but only a flattening. 



for the reasons discussed below. 

As found by many researchers (Waxman 1997, Panaitescu & 
Meszaros 1998, Sari 1998, Granot, Piran, & Sai'i 1999), at high 
frequencies, a spherical, relativistic fireball appears as a disk 
whose brightness surface increases toward the edge. An in- 
creasing magnification through lensing requires that the source 
is initially outside the Einstein disk. The peak magnification is 
obtained when, due to the source expansion (eqs. |jl|] and [^]), 
its edge crosses the Einstein disk. A significant magnification 
for a short duration (1-2 days for 000301c) constrains the size 
of the Einstein disk in two ways: it cannot be much smaller than 
the source size at the time of maximum lensing, as in this case 
the magnification would be weak, but it cannot be larger than 
it either, because the lensing event would last longer than the 
time when the maximum magnification is produced (3^ days 
for 00301c). For the jet models with microlensing presented in 
Figure 2^, the evolution of the source size and of the surface 
brightness distribution are such that the best fits are obtained 
if, at the time of maximum magnification 3.5 days), the af- 
terglow apparent size is twice larger than the Einstein disk. For 
this source-lens geometry, the source is too large for microlens- 
ing to yield a peaked light-curve. 

Garnavich et al. (2000b) have obtained peaked lensed light- 
curves by approximating the source image with an annulus, 
whose fractional width (16%) was obtained by fitting the data. 
Granot & Loeb (2001) have calculated microlensed light-curves 
using the surface brightness distribution of a spherical, relativis- 
tic afterglow, and have found a peak magnification factor ^ 2, 
supporting the microlensing model of Garnavich et al. (2000) 
for 000301c. We note, however, that in 000301c, the optical 
brightening occurs at times when jet effects are important, thus 
the source geometry and dynamics are not those expected for a 
spherical fireball. According to loka & Nakamura (2001), the 
surface brightness distribution of a spreading jet is rather uni- 
form. We also note that, at the time of optical brightening, the 
jets of Figure 3 and 4 are not very relativistic (F ^ several), 
which makes the surface brightness distribution even more uni- 
form. Therefore, in our modeling of 000301c, jet effects and 
the mildly relativistic motion of the source lead to an afterglow 
whose image at optical frequencies is not like the narrow ring 
envisaged by Garnavich et al. (2000b), but more like a disk. 

5. CONCLUSIONS 

The analysis presented in §^ shows that the "traditional" jet 
model cannot reconcile the late time decay indices of the 8.5 
GHz and i?-band light-curves of 000301c, and cannot account 
for the magnitude of the steepening exhibited by decay of its 
optical emission. We attribute the latter to the passage of a 
spectral break ly^ associated with a steepening of the power-law 
electron distribution at high Lorentz factors. This also provides 
a natural explanation for the curvature of the optical spectrum 
observed by Jensen et al. (2001) without the help of host extinc- 
tion, and for the possible softening of the NIR-optical spectrum 
implied by the results of Jensen et al. (2001) and Rhoads & 
Fruchter (2001), at about the same time when the steepening of 
the optical fall-off of 000301c occurred. Note that a chromatic 
evolution of the afterglow emission over less than a decade in 
frequency cannot be due to microlensing or to jet effects. 

A break in the electron distribution was also used by Li 
& Chevalier (2001) to model the emission of 000301c in the 
framework of spherical ejecta interacting with a wind-like ex- 
ternal medium, and by Panaitescu & Kumar (2001) for the af- 
terglow of GRB 991216, where its was found that a jet model 
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with a broken power-law injected electron distribution accom- 
modates the data better than the simple jet model. Although the 
optical light-curve breaks shown in Figures 3 and 4 are mostly 
due to the i^^-passage, the collimation of ejecta is still present 
in the model, to obtain consistency between the hard electron 
index p ^1.5 and the t^^ late radio decay. 

We find marginally acceptable fits only for a jet interacting 
with a homogeneous medium, and after excluding the 3.0^.3 
days data, as we could not model the brightening of 000301c 
during this time using microlensing. From the variation of 
around its minimum, we arrive at the following 95% confidence 
level intervals for a single parameter: 

Eo = (1.8 - 3.6) X lO^" erg , 6*0 = 10.8 - 14.6 deg , 

n=13-38cm-3, £3 = 0.032 - 0.15, 

£e = 0.030- 0.065 , p= 1.41 -1.59, 

6 = 0.34- 0.54, (5p = 0.9-1.3. 

These results are given only to show the statistical uncertainties 
of the model parameters, and should not be considered abso- 
lute limits. For such limits one would need to know the "model 
error bars", which may be larger than the above uncertainties. 

Consistency between the expected decay of jet light-curves 
and that seen in 000301c at 8.5 GHz after 30 days requires 
a hard electron distribution, with p ^ 1.6, similar to the in- 
dices obtained by MaUcov (1999) for Fermi acceleration in the 
limit when particles acquire a significant fraction of the shock 



energy, and by Ellison, Jones, & Reynolds (1990) and Baring 
(2000) for particle diffusion driven by large angle scatterings. 

We note that the time when the afterglow steepening is seen 
constrains the break frequency and determines the electron 
fractional energy e. In light of this, is interesting to note that for 
000301c we find an electron fractional energy close to equipar- 
tition, providing thus a natural reason for steepening of the elec- 
tron distribution at high energies. 

The model lensed light-curves we obtain show that mi- 
crolensing by can yield a significant magnification but does 
not provide a good fit to the observations, failing to produce 
a peak at maximum lensing. The reason is that the afterglow 
surface brightness is not sufficiently concentrated toward the 
edge of the visible disk to yield a light-curve peak when this 
edge crosses the Einstein disk of the lens. We emphasize that 
this conclusion is subject to the correctness of the assumptions 
made in our treatment of the jet model. 

Evidently, we do not rule out microlensing of 000301c within 
the framework of other afterglow models (e.g. non-uniform 
jets). An inhomogeneity in the external medium or a pile-up 
of electrons (Protheroe & Stanev 1999) at the 7* break, above 
which electrons are inefficiently accelerated, could in principle 
explain the brightening that preceded the steepening of the op- 
tical emission of 000301c. 



It is a pleasure to thank Pawan Kumar (IAS) for useful sug- 
gestions and comments on this work 
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Fig. 1 . — Spectrum of the 000301c afterglow att = 5.0 day constructed from the measurements closest to this epoch and linear inteipolations at low frequencies or 
using the fit of Rhoads & Fruchter (2001) in the NIR-optical. Optical fluxes are corrected for Galactic extinction of E{B — V) = 0.053. Downward triangles indicate 
1(T and 2cr upper limits (note that Berger et al. 2000 report a flux of 2.85 ± 0.95 mJy at 100 GHz an 4.26 days; the upper limits shown here are their measurement at 
5.09 days). The low flux at 10^^ Hz, representing the UV measurement of Smette et al. (2001), could be due to intergalactic Lya absorption. The solid, dashed and 
dotted fines represents the spectra obtained with the jet models described in Figures 2, 3, and 4, respectively. 
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Fig. 2. — Best fit models for the 000301c afterglow with a soft electron distrihution, homogeneous external medium and cooling frequency above optical. 
Triangles show la and 2a upper limits. Dotted hnes: jet model with parameters (j |2. l| ) Eq = 1.9 X lO""^^ erg, 9o = 3.1°, n = 0.015 cm~'^, ee = 0.050, 
eb = 4.1 X 10"*, p = 2.56, lens parameters Oj^s = 4.5 ^tas and 9e = 2.7/ias, and = 680 for 132 df. Dashed Hnes: jet model without microlensing, 
parameters Eq = 1.8 X 10^^ erg, 6»o = 2.8°, n = 0.011 crQ-^, ee = 0.047, £s = 4.4 X 10""', p = 2.53, and = 480 for 98 df, excluding the 3.0^.3 day 
data. Vertical dotted segments illustrate the modulation amplitude due to interstellar scattering. Optical light-curves have been displaced vertically by the indicated 
factors. The afterglow is at redshift z = 2.03 (Castro et al. 2000, Smette et al. 2001). We assumed Hq = 65 kms^^Mpc^i, go = 0.1, and A = 0. 
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Fig. 3. — Best fit models with a hard electron distribution, cooling frequency below optical, and /lomogeneous external medium. The jet model shown with dotted 
lines includes microlensing, all data, and has = 225 for 130 df, Eq = 2.0 X 10'™ erg, 6»o = 13.6°, n = 15 cm"^, es = 0.18, Ee = 0.045, p = 1.33, 
e = 0.33, Sp = 1.1, Bji^s = 0.90 ^as and 8^ = 0.38 fj,as. Dashed lines are for a jet model without microlensing and the 3.0—4.3 day measurements, and has 
= 120 for 96 df, Eo = 2.7 x 10'™ erg, 6»o = 11.1°, n = 17 cm-^, eb = 0.072, Se = 0.037, p = 1.49, e = 0.31, 5p = 1.1. The vertical lines show the 
interstellar scintillation amplitude. 
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Fig. 4. — Best fit models with a hard electron distribution, cooling frequency below optical, and wind-lilie external medium. Dotted lines represent the model with 
microlensing, all data, x^ = 330 for 130 df, and parameters Eq = 1.5 X 10^0 erg, = 7.8°, A* = 0.94, eb = 0.25, Ee = 0.022, p = 1.36, e = 0.32, 
Sp = 1.6, Ols = 0.62 /las and Be = 0.29 fias. Dashed lines are for a model without microlensing, excluding the 3.0-4.3 day data, and has x^ = 140 for 96 df, 
Eo = 1.9 X 10^0 erg, 6»o = 10.6°, A, = 0.90, eb = 0.22, Ee = 0.047, p = 1.36, e = 0.47, Sp = 1.7. 



